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Abstract 

Singularities in the dark energy universe are dis- 
cussed, assuming that there is a bulk viscosity in the 
cosmic fluid. In particular, it is shown how the physi- 
cally natural assumption of letting the bulk viscosity be 
proportional to the scalar expansion in a spatially flat 
FRW universe can drive the fluid into the phantom re- 
gion (w < —1), even if lies in the quintessence region 
(w > —1) in the non-viscous case. 
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1 Introduction 

Recent astrophysical data indicate the presence of a mysterious 
kind of energy (ideal fluid with negative pressure), contributing 
about 70% of the total energy of the universe. It is quite possible 
that the equation of state parameter w for dark energy is less 
than -1. If this is so, the universe shows some very strange 
properties such as the future finite singularity called Big Rip pQ. 
In turn, this can lead to the occurrence of negative entropy 
One may expect that natural effects, most likely of a quantum 
mechanical origin, may prevent the Big Rip, as was shown in 
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Ref. j3J. Dependent on which sorts of dark energy are present, 
various types of singularities appear; for a classification, see 
Ref. 0J. In the present letter, we will consider the role of a 
bulk viscosity ( in a universe having a Big Rip singularity We 
take the universe to be spatially flat, with cosmo logical constant 
equal to zero. The shear viscosity 77 will be put equal to zero, 
in conformity with usual practice. The latter assumption is 
however a non-trivial point, since the shear viscosity is usually 
so much greater than the bulk viscosity. At least, this is so in 
the early universe. For instance, in the plasma era after the time 
of recombination (T ~ 4000 K), the ratio r]/( as calculated from 
kinetic theory is as large as about 10 12 jH]. Therefore, even a 
slight anisotropy in the cosmic fluid would easily outweigh the 
influence from the bulk viscosity. We will below be concerned 
with the late universe and not the early one. However, we ought 
to bear in mind that the neglect of rj most likely rests on the 
tacit assumption about an extreme fine tuning of the spatial 
isotropy of the fluid. 

In the next section we present the formalism for the FRW 
universe when the equation of state is of a general form, w = 
w(p). Similarly, we admit a general form for the bulk viscosity, 
C — C(p)- Thereafter we discuss some special cases. Of main 
interest appears to be the case where ( is proportional to the 
scalar expansion 9, ( = r 9, r being a constant. Then, it turns 
out that the barrier w = — 1 between the quintessence region 
(w > —1) and the phantom region (w < —1) can be crossed, 
as a consequence of the bulk viscosity. To our knowledge, this 
has not been pointed out before. The case ( = r 6, implying 
proportionality of the bulk viscosity to the divergence of the 
fluid's velocity vector, is physical natural, and has been consid- 
ered earlier in an astrophysical context. Cf., for instance, the 
review article of Gr0n |B]. 

2 General formalism 

We assume the standard FRW metric, 

ds 2 = -dt 2 + a 2 {t) dx 2 , (1) 
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and put the spatial curvature k as well as the cosmological con- 
stant A equal to zero. The Hubble parameter is H = a/ a, and 
the scalar expansion is 9 = U^ tll = 3H, being the four- 
velocity of the fluid. We assume the equation of state in the 
form 

P = w (p)p> ( 2 ) 

where the thermodynamical variable w(p) is an arbitrary func- 
tion of the density. Similarly, we assume that the bulk viscosity 
is arbitrary, ( = ((p). The Friedmann equations are (cf., for 
instance, Ref. [5]) 

6 2 = 24vrGp, , (3) 
p + (p + p)9 = (9 2 , (4) 

-+ 6 ^ = -^G(p-(6). (5) 
a io 

The effective pressure is p = p — (9. On thermodynamical 
grounds, in conventional physics ( has to be a positive quantity. 
This is a consequence of the positive sign of the entropy change 
in an irreversible process (cf., for instance, Ref. [7]). We shall 
assume that ( > also here, although although strictly speaking 
the usual thermodynamical relationships ought to be taken with 
some care when dealing with bizarre systems like a phantom 
fluid. The effect of the bulk viscosity is accordingly to reduce the 
thermodynamic pressure in the fluid. (The viscous Friedmann 
equations above were considered in a Cardy-Verlinde entropy 
context in Ref. jH].) Note that due to the time dependence of the 
viscosity the Friedmann equations above may also be interpreted 
as coming from a modification of gravity (for a recent discussion 
and a list of references, see Refs. |Hj and [TUj). 

From the above equations we derive the following differential 
equation for the scalar expansion: 

9 + hl + w)9 2 - 12tiG(9 = (6) 

•Li 

which, in view of the relationship 9 = y&nGp/ y/p, can be writ- 
ten as a differential equation for p. Since we will be interested 
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in the region around w = — 1, it is convenient, following Nojiri 
et al. jl], to introduce the function f(p) defined by 



l+w(p) = -f(p)/p, (7) 

implying that f(p) = for a "vacuum" fluid. Then, the equation 
takes the form 

p - y/2^Gp~f(p) - 24vrGC(p)p = 0. (8) 

We shall be interested in the development of the late universe, 
from t = to onwards. For simplicity, we put to = 0. The corre- 
sponding starting value of p will be denoted by p . From Eq. (jHJ) 
we obtain 



t 



dp 0) 



This is the general relation between the cosmological time t (> 
0) and the density p. 



3 The case when w is a constant 

This is thermodynamically the simplest case. Let us put 

f(p) = up, (io) 

where a is a constant. This means that p = wp = —(1 + a) p. 
We investigate in the following some different choices for (. 
(i) C — 0. In this non- viscous case Eq. yields 

(= i l(-L-± 

V247TG a Vv/Po y/p 

Thus, if a > a finite value of the time t is compatible with a 
Big Rip singularity (p = oo). This is the conventional phantom 
case, corresponding to io < —1. The variation of the density 
with time is 2 

pit) = p (i - \ote Q t\ , (12) 
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where the initial scalar expansion is 



9 = v^vrGpo 
Similarly, the scalar expansion is 

1 



9{t) 

and the scale factor is 



1 a9nt 

2 



a(t) — a ( 1 — -a9 t 



-2/3o 



(13) 



(14) 



(15) 



One may note here that pj po = (a/ao) 3a . Thus all the quantities 
pit), 9{t), a(t) diverge at the Big Rip. The above expressions are 
in agreement with earlier results PQE|. It ought to be mentioned 
that for w < — 1 the dark energy behaves in close analogy with 
quantum fields [TTj . 

(ii) ( equal to a constant. Equation © yields in this case 



dp 



po p 



3/2 



1 + v / 247GC/(a v / P) 



(16) 



It is however in this case most convenient to go back to Eq. 
and solve it with respect to 9 (cf. Ref. jHj): 



9{t) 

where 

Correspondingly, 

a(t) = ao 
From Eq. © then 

P(t) = Po 



t/t c 



1 - \a9 Q t c {e^ - 1) : 
t c = (UttGC)- 1 . 



1 - -a9 t c (e t/tc - 1) 



(17) 



-2/3a 



e 2t/t c 

[1 - \a6 Q t c {e^ - 1)] 
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(19) 



(20) 



Thus, at the time t = t s , where 



t x = t c In 



1 + 



a0 t c 



(21) 



there occurs a Big Rip; both 6{t),a{t) and p{t) diverge. As 
before, a has to be a positive quantity, corresponding to w = 
— 1 — a < — 1. When ( — > 0, the results of the previous case (i) 
are recovered. 

(in) The case ( = t9. This is the most interesting case. As 
C is assumed to be positive, the constant r has to be positive. 
In view of Eq. (jSJ) we can alternatively write the bulk viscosity 
function of the density: 



C(p) = T^/24nGp. (22) 

Equation Q yields in this case 

1 2 / 1 1 \ 

t = -== — — — ^ (—- — ). (23) 



From this condition it follows that the condition for a Big Rip 
(p = oo) to occur in a finite time t is that the prefactor is 
positive, 

a + 2AttGt > 0. (24) 

This is the most important result of the present paper. Even 
if we start from a situation where w > — 1 (i.e., a < 0), corre- 
sponding to the quintessence region for an ideal fluid, the pres- 
ence of a sufficiently large bulk viscosity will make the condition 
(J2IJ) satisfied and thus drive the fluid into the Big Rip singular- 
ity. Somewhat surprisingly, it is physically the reduction of the 
thermodynamical pressure generated by the viscosity which in 
turn causes the barrier w = — 1 to be crossed. 

The system is actually most easily analyzed by going back 
to Eq. (jHJ) for the scalar expansion. Introducing the effective 
new parameter 

d = a + 24ttGt, (25) 
we can write the equation as 

0-}-a9 2 = O. (26) 
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From this it follows that p(t), 0(t) and a(t) can actually be found 
from the non- viscous expressions (|12j ) - (jl5p . only with the re- 
placement a — > a. 

In general, the existence of viscosity coefficients in a fluid is 
due to the thermodynamic irreversibility of the motion. If the 
deviation from reversibility is small, the momentum transfer 
between various parts of the fluid (i.e., the stress tensor) can 
be taken to be linearly dependent on the velocity derivatives 
dkUi. This is the conventional case, corresponding to constant 
viscosity coefficients. The present case ( = t9 means that we 
go step further; it means that the momentum transfer involves 
second order quantities in the deviation from reversibility, still 
maintaining the scalar property of Q. As mentioned above, the 
ansatz ( = t8 has been considered earlier in a cosmological 
context 0. 

4 Remarks on the case w = w(p) 

In this general case we have to go back to the expression Q. 
As the ansatz ( = t9 appears to be that of main interest, we 
consider this case first. From Eq. Q we get 



Thus, if f(p) — > ap in the limit of large p, the Big Rip singularity 
is allowed. The necessary condition on the values of a and r are 
again given by Eq. (j2U). 

Another choice for f(p) that may seem natural, is to take 



for all p, where A and (3 are constants (cf. Refs. jU H21 EH])- 
Thus p = — p — ApP . We shall here take A and (3 to be positive. 
From Eq. (f2*7jl we get 




(27) 



Hp) = Ap- 



(28) 




(29) 
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If (3 > 1 the influence from viscosity is seen to fade away for 
large p, and the Big Rip is allowed. If (3 < 1 the viscosity term 
becomes dominant for large p, and if we perform the integration 
over a region p G [po, p] where p is large, we get 

2 1/1 1 \ 

f W (24ttG)3/2 7 ^ " ■ (30) 

We can thus have Big Rip also when j3 < 1. Close to the 
singularity, the time t depends on viscosity through the factor 
r but is independent of A and (3. Recall that the formalism so 
far in this section assumes that £ = t9. 

The final point that we shall deal with, is the following ques- 
tion: what form for the function f(p) does the ansatz 

fl W = a ° (vtt) (31) 

correspond to? Here n is a positive constant and t s is the Big Rip 
time. The ansatz (|H1)I is simple and suggestive, and has been 
dealt with before in the non- viscous case [I] . From Eqs. (J3T|) 
and © we get 

3n 2 / 1 1 \ 2 
PW = ^ 7 + ^- , (33) 



8ttG V t t. 



leading to 




m = ± 2 P ^JL-^l^ P> . (34) 

The doubly-valued equation of state is typical for a first-order 
transition (cf. a more extensive discussion on this point in 
Ref. 03). These equations are the same as in the non- viscous 
case. But there is a difference arising from Eq. (@J): we get 

1/2 
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showing the influence from viscosity in the last term. Note that 
C = C(p) is nere an arbitrary function. Thus, the initial simple 
ansatz (|3Tj) for a(t) leads to a rather complicated form for f(p). 

It would be interesting to generalize the present work tak- 
ing into account quantum effects together with viscosity. This 
may be done in close analogy with previous attempts in this 
direction, for an asymptotically de Sitter universe j!4j . 
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